Soliton topology versus discrete symmetry in optical lattices 
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We address the existence of vortex solitons supported by azimuthally modulated lattices and reveal 
how the global lattice discrete symmetry has fundamental implications on the possible topological 
charges of solitons. We set a general "charge rule" using group-theory techniques, which holds for 
all lattices belonging to a given symmetry group. Focusing in the case of Bessel lattices allows us 
to derive also a overall stability rule for the allowed vortex solitons. 



Vortex solitons in nonlinear systems (for a review see 
0) characterized by a discrete symmetry have been nu- 
merically predicted in two-dimensional arrays of evanes- 
cently coupled waveguides 0, harmonic refractive index 
gratings imprinted in cubic media 0, as well as in pho- 
tonic crystal fibers with defects 0. Because of the im- 
printed refractive-index modulation, such vortices can be 
made stable in contrast to their ring-shaped counterparts 
in uniform focusing media. Recently, vortices having unit 
topological charge have been experimentally observed in 
optically-induced lattices in photorefractive media ||. 
The very refractive index modulation causing the stabi- 
lization of vortex solitons simultaneously imposes restric- 
tion on the possible topological charges of the vortices 
dictated by the finite order of allowed discrete rotations 
. A corollary of such result is that the maximum charge 
of stable symmetric vortex in two-dimensional square lat- 
tices is one. 

However, square lattices are just one particular exam- 
ple of guiding structures accessible for experimental ex- 
ploration. Another interesting class of such structures 
with a new global rotational symmetry is constituted by 
azimuthally modulated lattices, also offering a wealth of 
new opportunities. For example, in such lattices the or- 
der of rotational symmetry may be higher than 4, in con- 
trast to square lattices, a property that has direct impli- 
cations in the possible topological charges of symmetric 
vortex solitons supported by such lattices. 

In this Letter we explore the connection existing be- 
tween the lattice discrete symmetry and the topology of 
the allowed vortex solitons by means of a general group- 
theory approach. We find that azimuthally modulated 
lattices imprinted in focusing medium can support sym- 
metric vortex solitons carrying phase dislocations with 
topological indices higher that one. The higher the sym- 
metry order of the lattice, the higher the allowed vortex 
topological charge. We also address the stability of the al- 
lowed vortex soliton families, taking as a particular exam- 
ple the case of azimuthally modulated optically-induced 
Bessel lattices. 

Optical lattice induction in anisotropic nonlinear ma- 
terials introduced in Q opens broad prospects for cre- 



ation of reconfigurable refractive index landscapes with 
different types of nondiffracting beams, including Bessel 
beams || . Accurate approximations of Bessel beams can 
be generated experimentally in a number of ways. Known 
techniques include illumination of annular slit in the focal 
plane of a lens, conical axicons, as well as more compli- 
cated interferometric and holographic techniques (see 
and references therein). Thus, azimuthally modulated 
lattices of any desired order can be optically induced by 
higher-order Bessel beams 0. Because of the diffrac- 
tionless nature of Bessel beams, they are to be launched 
collincarly along the anisotropic nonlinear material (e.g., 
photorefractive crystal) with a polarization orthogonal 
to the soliton beams, to make use of incoherent vecto- 
rial interactions Q. The concept can be extended to all 
relevant physical settings, including Bose-Einstein con- 
densates. 

Our starting point is the paraxial nonlinear equation 
for the complex field amplitude q describing the propa- 
gation of light in azimuthally modulated lattice: 



drj 2 



q-q\q\ 2 -pR{ V ,Qq. (1) 



The longitudinal £ and transverse 77, ( coordinates are 
scaled to the diffraction length and to the input beam 
width, respectively. The parameter p accounts for the 
depth of the refractive index profile, whereas the func- 
tion R(r), C) describes the lattice profile. In the par- 
ticular case of optical lattice induced by higher-order 
Bessel beam the refractive index profile features beam's 
intensity distribution Rn(v> = ■^n[(2&iin) 1/ ' ] cos 2 (n^), 
where r = (rj 2 + C 2 ) 1 ^ 4> l& the azimuth angle, and &i; n 
is the parameter that sets the transverse lattice scale. In 
the particular case of optical lattice induction in SBN 
crystal biased with dc electric field ~ 10 5 V/m, for 
laser beams with width 10 /im the propagation distance 
£ ~ 1 corresponds to I mm of actual crystal length, while 
amplitude q ~ 1 corresponds to peak intensity about 
50mW/cm 2 . The intensity of lattice-creating wave is of 
the same order and can be varied to tune p in a broad 
range. 

The realization of rotational symmetry in azimuthally 



modulated lattices is different depending on the order 
of modulation n. Lattices of lowest-order with n = 
show perfect continuous rotational symmetry since the 
function R (r],() = Ro( r ) nas no dependence on <j> 
In terms of group theory, the symmetry group of a lat- 
tice of zero order is 0(2). The azimuthal modulation 
R n (r, 4>) ~ cos 2 (ruf)) changes the symmetry group asso- 
ciated to the rotation transformations. The rotational 
symmetry group of a lattice of order n is given by the 
discrete point-symmetry group C-x n ,v, corresponding to 
discrete rotations of angle e n — ixjn with respect to a 
rotation axis perpendicular to the plane and intersecting 
it at the origin (R n (r,4> + e n ) — R n (r,(/))) as well as to 
specular reflections with respect to a number of planes 
containing the rotation axis [loj . This fact has strong 
implications in the form of possible symmetric vortex so- 
lutions of Eq.(IJ q(r], C,£) = [u(r], () + iv(rj, £)] exp(ife£), 
where u and v are real and imaginary parts, respectively, 
& is a propagation constant, and m describes topologi- 
cal winding number of complex field q that can be de- 
fined by the circulation of the gradient of the field phase 
arctan(v / u) around the singularity at 77, £ = 0. 

Individual vortices are characterized by a phase singu- 
larity, accompanied by a single point of zero amplitude. 
Discrete rotations are defined with respect to this point, 
so that all our group theory arguments will apply to in- 
dividual vortices and not to bound states or clusters of 
individual vortices owning intrincated phases featuring 
edges and multiple phase singularities. In this context, 
individual vortex solutions appear as doubly-degenerated 
pairs belonging to the two-dimensional representations of 
C2n.11 (or of its subgroups) and they are characterized by 
the index representation m. The charge of these solu- 
tions is exactly provided by the index representation m 
. For every value of m one finds a degenerated vortex- 
antivortex pair, with charges m and — m, respectively. 
Since m is limited by symmetry constraints, an upper 
bound for the values of permitted vortex charges is es- 
tablished if the system enjoys a rotational symmetry of 
finite order N @: < m < N/2 (for even N). Since 
the symmetry group of a lattice of order n is Ci n ,vi its 
symmetry order is N — 2n and thus even. Consequently, 
we obtain one of the central results of this Letter in the 
form of a "charge rule" for the allowed values of vortex 
charges in a lattice of order n: 



< m < n — 1 



(2) 



We stress that this rule is obtained on the basis of gen- 
eral symmetry arguments and is applicable for a whole 
class of azimuthally modulated lattices, irrespectively of 
the details of their local shape and thus method of re- 
alization. Some implications of [3 are readily apparent. 
For example, azimuthally modulated lattices of high or- 
ders allow generation of symmetric vortex solitons with 
topological charges higher than one, in contrast to square 
lattices. 




Figure 1: Amplitude (left column) and phase (right column) 
distributions for the vortex solitons with charges 1 (a), 2 (b), 
and 3 (c), supported by the fourth-order Bessel lattice at 
p = 28. All solitons correspond to b — 1.6. In left column 
bright regions correspond to high light intensities, dark re- 
gions correspond to low intensities. 



We performed a comprehensive numerical analysis of 
vortex soliton solutions of Eq. Q for the case of optically 
induced Bessel lattices to confirm the "charge rule". We 
searched for symmetric stationary vortex profiles with a 
relaxation method using as initial guess functions with 
a phase profile exp(im<j)) carrying phase-singularity with 
charge m. The numerical results confirmed the "charge 
rule" in all cases. A summary of this analysis is pre- 
sented in Table Q] No higher-order vortex solitons have 
been found to exist above the maximum limit given by 
Eq. (J2J. The amplitude and phase distributions for al- 
lowed vortices with charges 1, 2 and 3 in a Bessel lattice of 
fourth order are shown in Fig. ^ The presence of focus- 
ing nonlinearity tends to increment the localization effect 
generating the typical patterns of bright well-localized 
spots. Notice, that in sufficiently deep lattices at fixed b 
and p the bright spots in the vortex intensity distribution 
become more pronounced with increase of vortex charge 
m, while radius of the vortex remains almost unchanged 
since it is mainly determined by transverse lattice scale. 
This is in contrast to vortices in uniform medium that 
broaden notably with increase of m at fixed b. 

The properties of vortex solitons supported by fourth- 
order azimuthally modulated Bessel lattice are summa- 
rized in Figs. |2£a)-(d). The vortex soliton is character- 
ized by its energy flow U = J j \q\ 2 drjd(. The depen- 
dence U(b) for a lowest-order vortex soliton with m = 1 
is shown in Fig. |2Ja). The different behavior of the U(b) 
curve for a shallow (p = 4) or deep (p = 14) lattice in- 
dicates the importance of the confining properties of the 
Bessel lattice. The energy flow is a non-monotonic func- 
tion of the propagation constant for small lattice depths 
(p = 4) whereas for large values of p (approximately 
more or equal than 6) it monotonically increases with 
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Figure 2: (a) Energy flow U versus propagation constant b 
for a vortex soliton with m = 1. Propagation constant cut- 
off 6 co versus lattice depth p for vortex solitons with m — 1 
(b) and m = 3 (c). (d) Real part of perturbation growth 
rate 8 versus propagation constant at p — 28. Panels (a)-(d) 
correspond to fourth-order Bessel lattice. Real part of per- 
turbation growth rate versus propagation constant for vortex 
solitons in (e) third- and (f) fifth-order lattices at p = 28. 

b. The deeper is the lattice the more pronounced is the 
azimuthal modulation of the vortex intensity profile. 

In contrast, to = 1 vortices in shallow lattices may un- 
dertake an abrupt change of behavior at a critical value 
of the propagation constant (see the inflection point of 
the U(b) curve at p = 4). With increase of energy flow , 
bright spots forming such vortices may eventually merge 
into a single weakly modulated bright ring, since the con- 
fining effect of the shallow Bessel lattice is not sufficient 
to ensure the trapping of light in the higher refractive 
index regions. At a given value of p the energy flow of 
the vortex soliton vanishes at certain cutoff b co on prop- 
agation constant (Fig. H^a)). Close to the cutoff, vortex 
with charge to = 1 broadens drastically in shallow lat- 
tices, while in deep lattices brigth well-localized spots are 
always resolvable in vortex intensity profile. As shown in 
Fig. Elb), for vortex with charge to = 1 the cutoff is 
a monotonically increasing function of the lattice depth 
p. For higher-order vortices with charges m > 1 the de- 
pendence b co (p) may be discontinuous (Fig. [2Jc)). Thus, 
in shallow lattices higher-order vortices cease to exist in 
the cut-off without any topological transformation (left 
branch of the b co {p) curve), while in deep enough lat- 
tices they drastically broaden in the cutoff where vortex 
energy flow vanishes (right branch). 

Next we conducted a detailed linear stability analysis 
for the allowed vortex soliton families. We searched for 



perturbed solutions of Eq.QJ in the form q(r],(:0 = 
[u(r), C) + u p (rj, C, + iv(r), Q + iv p (r), (, £)] exp(i6£), 
where u p and v p are real and imaginary parts of pertur- 
bation, respectively. The linearized evolution equations 
for perturbation components are 

du 1 

~ajr = — -^A±v p + bvp - [2uvu p + (3v 2 + u 2 )v p ] - pRv p 



— = --A±Up + bu p - \2uvVp + (3« 2 + v 2 )u p ] —pRu p 

where Aj_ stands for the transverse Laplacian. We solved 
this system numerically by means of a split-step Fourier 
method with noisy initial conditions to get perturbation 
profiles and their growth rates. 

The results of the stability analysis are summarized in 
Figs. Efd)-(f) and in Table [I] for lattices with n up to 
6. Thus, in fourth- and fifth-order Bessel lattices (Figs. 
I3d) and|2ff)), vortices with m = 1 and 2 were found to 
be unstable in the entire domain of their existence, while 
in third-order lattice the unstable vortex carries charge 
rn = 1 (Fig. Efe)). Both exponential and oscillatory in- 
stabilities are encountered for such vortices. In contrast, 
third-order lattice can support stable vortex with charge 
m = 2, fourth-order lattice supports stable vortex with 
m = 3, and fifth-order lattice supports stable vortices 
with to = 3 and 4. In shallow lattices the domains of 
existence of such vortices feature multiple stability and 
instability windows, while for deep enough lattices such 
vortices become stable in the entire domain of their ex- 
istence. The physical origin of such stabilization is the 
higher confinement of radiation in regions with higher 
refractive index that reduces the strength of nonlinear 
interactions of bright spots forming the vortex. In exper- 
iment vortex stabilization would become apparent upon 
gradual increase of intensity of the lattice-creating beam. 

Similar comprehensive linear stability analysis con- 
ducted for all lattices with orders n up to 20 with various 
depths < p < 100 enabled us to derive the important 
"stability rule" for vortex solitons, which may be stable 
(i.e., a necessary but not sufficient condition) only if the 
vortex topological charge satisfies the condition 

Ti 

- < to < n - 1 (3) 

with the exception for n = 2, when the only existing sym- 
metric vortex with charge to = 1 may be stable. Similarly 
to the cases of third-, fourth-, and fifth-order Bessel lat- 
tices (Figs. Eld)-(f)), vortices whose charges satisfy the 
condition |3 become stable in the entire domain of their 
existence for deep enough lattices. 

Notice that the phase variation S<fi between neighboring 
lattice maxima for vortices with charges verifies tt/2 < 
8(j) < 7r. On intuitive grounds, this difference is consistent 
with a repulsive interaction between bright spots forming 
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Lattice order 


Available charges and stability status 


n = 2 


m = 1 
stable 










n = 3 


m = 1 
unstable 


m = 2 
stable 








n = 4 


m = 1 
unstable 


m = 2 
unstable 


m = 3 
stable 






n = 5 


m = 1 
unstable 


m = 2 
unstable 


m = 3 
stable 


m = 4 
stable 




n = 6 


m = 1 
unstable 


m = 2 
unstable 


m = 3 
unstable 


m = 4 
stable 


m = 5 
stable 



Table I: Table showing the available charges and the stability 
status of vortex solutions for different lattice orders. The 
stable status means that it is possible to stabilize the vortex 
under consideration by a suitable election of the Bessel lattice 
parameters. 



;=0 

ft 
• # 




£=20 
• 


•• 


• 

• 





* ■ • * 



Figure 3: Propagation dynamics of vortex solitons with m = 3 
(a) and m — 1 (b) supported by fourth-order Bessel lattice in 



the presence of white input noise with variance a no 



0.01. 



Bright regions correspond to high light intensities, dark re- 
gions correspond to low intensities. Both solitons correspond 
to b = 1.6 and p = 28. 



the vortex (Fig. [IJ, that are compensated by the lattice 
and lead to stable vortex propagation, similarly to the 
case of multipole beams ||. 

Direct numerical integration of Eq. Q with input con- 
ditions g(ry,C,C= 0) = w(r], Q[l + p(rj, C)], where io(»7,C) 
is the stationary solution and p(rj, £) is the white noise 
with variance cr„ oise = 0.01, fully confirmed results of 
linear stability analysis. In the presence of noisy per- 
turbations, vortices whose charges satisfy the "stability 
rule" propagate undistorted over hundreds of diffraction 
lengths, while unstable representatives of vortex soliton 
families are rapidly destroyed upon propagation (FigJHJ. 
The decay of unstable vortex is accompanied by progres- 
sive increase of intensity oscillations in neighboring bright 
spots, until only several spots remain in the output pat- 
tern. 



In conclusion, we set a general connection between the 
discrete rotational symmetry of azymuthally modulated 
lattices and the topological winding number of the al- 
lowed symmetric vortex solitons. We derived a "charge 
rule" and "stability rule"; the "charge rule" is intended 
to be general for lattices with a given global symmetry 
order and does not depend on the way in which lattice 
is created (optical induction in anisotropic nonlinear ma- 
terial or in Bose-Einstein condensates, direct fabrication 
in the case of photonic crystal, or other methods), while 
the details of stability of the allowed vortex families are 
expected to depend on the local properties of the partic- 
ular lattice considered. The rules predict, for example, 
that in contrast to square lattices, suitable azimuthally 
modulated lattices support stable symmetric vortex soli- 
tons with higher-order topological charges. Because of 
the general nature of the geometrical method used in 
our derivation, we anticipate that a similar charge rule 
should hold for scalar soliton systems with comparable 
azymuthal symmetry in a variety of physical settings. 
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